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In this paper we investigate the behaviour of a fully quantum mechanical system consisting of a 
mesoscopic SQUID ring coupled to one or two electromagnetic field modes. We show that we can use 
a static magnetic flux threading the SQUID ring to control the transfer of energy, the entanglement 
and the statistical properties of the fields coupled to the ring. We also demonstrate that at, and 
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in the coupled system, and the energy transfer between its components, is strong. 
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I. INTRODUCTION 



In an earlier publication^ we considered the interaction 
of a quantum mechanical SQUID ring (a single Joseph- 
son weak link, capacitance C^, enclosed by a thick su- 
perconducting ring, inductance Ag) with a classical elec- 
tromagnetic (em) field. Using quasi-classical Floquet 
|.]^gpj.jj2j3j4j5j6j7 solve the time dependent Schrodinger 
equation (TDSE) for the SQUID ring, we were able to 
show that the ring-field interaction could be very highly 
non-pcrturbative in nature. In essence this is due to the 
ring Hamiltonian^ containing a cosine term (the Joseph- 
son coupling energy) which can generate non-linearities 
to all orders. In addition, this Hamiltonian and its so- 
lutions are $0 {— /i/2e)-periodic in the external static 
magnetic flux {^xstat) applied to the ring. This quan- 
tum non-linearity ensures that energy exchange between 
the field and the ring is dominated by multiphoton ab- 
sorption (and emission) processesi. As we have demon- 
strated, this is the case even at modest field amplitudes 
and at frequencies much less than the separation between 
the ring energy levels {~^h). In this work we showed that 
these energy exchanges occurred over very small regions 
in the bias flux ^xstat- The values in ^xstat at which 
these exchanges take place are determined by the ring 
energy level structure and the field frequency (ci;e/27r) 
and fiux amplitude ($e)- To be precise, it is in the ex- 
change regions that the energy expectation value appears 
to jump (for example, using a two level model) between 
the time-averaged energies of the ground and first excited 
states of the ring. Each transition (exchange) region cor- 
responds to the separation between the ring eigenenergies 
equalling nhwe, n integer, leading to multiphoton absorp- 
tion, or emission, between the ring and the field. It is 
in these regions that the non-linear nature of the ring 
Hamiltonian becomes manifest and where strong (and 



non-perturbative) time dependent superpositions occur 
between the original eigenstates of the ring. 

Currently there is a great deal of interest in using meso- 
scopic SQUID rings (and other weak link based circuits) 
in quantum technologies, for example, in quantum com- 
py^jj^g^jigjii^^ This interest has been stimulated by re- 
cent experimental work on probing quantum mechani- 
cal superposition states in Josephson weak link circuit 
system3i2*i^*iSii&, and even more so in the last year by 
reports of superposition states in SQUID ring9i£*iS*iS. 
It seems reasonable to assume that the theoretical de- 
scription of weak link systems interacting with em fields 
(classical and quantum mechanical) is likely to be of 
great importance in the development of any future su- 
perconducting quantum technologies. In this regard the 
very strong non-linear behaviour exhibited by a single 
weak link SQUID ring in the exchange regions, referred 
to above, may prove to be of great utility. In order to 
test this viewpoint we have recently considered, within 
a fully quantum mechanical framework, the interaction 
of a SQUID ring with an oscillator field modc^°, i.e. the 
simplest coupled system we could have chosen (see fig- 
ure^. We found that for the case of the em field in a 
coherent state the results derived from this quantum ap- 
proach compare very well with those obtained previously 
using quasi-classical Floquet theory. In both approaches 
the ring and the field mode only couple strongly together 
within the exchange regions, i.e over certain narrow re- 
gions in the bias flux ^xstat- This means that ^xstat can 
be used to control the coupling without losing superpo- 
sition coherence in the system. We note that this work 
relates to quantum optical interactions in few level atoms 
and to few level systems involving either (superconduct- 
ing) electron pairs or single electronsSiSiiSSiSiSi. 

These initial results for a two mode (ring -I- oscilla- 
tor) system have encouraged us to draw more parallels 
with quantum optics. Rather than simply consider the 
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FIG. 1: Block diagram of a SQUID ring coupled to a single 
em field mode of frequency uje (= 1/V CeAe) where the flux 
linkage factor, ring to field mode, is fies- Here, it is assumed 
that the temperature T is such that ksT <^ hjjjs,hu>e for a 
SQUID oscillator frequency Ws = 1/\/CIA7 • Also shown is a 
static bias magnetic flux ^xatat applied to the ring. 
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FIG. 2: Block diagram of a SQUID ring coupled to two 
em field modes of frequency ci^ei and tJe^ assuming ksT <^ 
hu)s, hujei, hu)e2, with flux linkage factors fieis and ^J.se2 be- 
tween, respectively, the first field mode and the ring and the 
ring and the second field mode; all else as for figure 



SQUID ring as an electronic device, we may also view it 
as a tunable, $o-periodic, non-linear medium to couple 
a system of quantum oscillators together. Regarded as 
a non-linear medium, there is a clear analogy to other 
non-linear quantum systems in the context of quantum 
optics. However, there are two crucial differences, both 
of which may be of great importance in future quantum 
technologies. First, unlike the SQUID ring, in quantum 
optical systems the medium usually displays a weak poly- 
nomial non-linearity, even in strong f[ci(f£25, 26,27^28,29^30^ 
Second, all of the properties of the SQUID (quantum or 
quasi-classical) are $o-periodic in bias flux. 

In this paper, our objective is to explore the con- 
sequences of the strong quantum non-linearity of the 
SQUID ring on the interaction, via the ring, of two os- 
cillator field modes. This arrangement is depicted in fig- 
ure 12 with the two held modes and the ring oscillator 
frequencies taken to be Uei/^TT, We2/27r and LOgjl-K = 
271-VAsC ' I'sspectively. As we shall see, the addition of 
the second field mode makes this a much more sophisti- 
cated and interesting system than the two mode system 
(ring-|-field oscillator) which was the subject of a recent 
publicationSSi, even though the computational demands 
that need to be met are very much greater. In this regard 
it is widely viewediSf^ that the SQUID ring, as a coher- 



ent quantum device, has many potential applications in 
the design, development and operation of quantum me- 
chanical circuits and quantum logic elements. In this pa- 
per we consider two aspects of the quantum behaviour of 
a SQUID ring which could have a serious impact in these 
areas, namely the transfer of entanglement and frequency 
conversion between em field modes via the quantum non- 
linearity of the ring. In this work we discuss frequency 
conversion and entanglement for just two field modes. 
However, if the non-linear aspects of quantum SQUID 
ring behaviour can be fully exploited more complicated 
operations could be envisaged. These may include using 
SQUID rings to couple/decouple entangled states in ex- 
tended qubit circuit structures and allow frequency con- 
version processes between field modes to be modulated, 
producing coherent pulse modulated signals. In our opin- 
ion, the combination of such strong non-linear proper- 
tics, coupled with $o-pcriodic external bias fiux control 
of this behaviour, makes the SQUID ring quite unique as 
a device for application in quantum technologies. Thus, 
although the following calculations are concerned with 
some of the basic consequences of the quantum interac- 
tion of em field modes with a SQUID ring, we also wish 
to emphasize the technological possibilities which may 
open up as these ring-field mode systems become more 
fully understood. 

In the work presented here we first consider briefly the 
two mode system, including a static bias flux ^xstat (fig- 
ureOJ . This allows us to relate the quasi-classical Floquet 
approach to the fully quantum mechanical treatment and 
demonstrate that our quantum model can produce con- 
sistent results. It also provides the background formal- 
ism for our main goal which is the study of two em field 
modes coupled through a SQUID ring. Since our purpose 
is to study the full quantum mechanics of the ring-field 
mode (1 or 2) system, we assume throughout that the 
operating temperature (T) is such that hwex.ei ^ ^bT 
hujs ^ ksT. This ensures that both the ring and field 
mode(s) behave quantum mechanically. We then con- 
sider the extended quantum circuit, the two oscillator 
field modes (A and £2) coupled through a SQUID ring 
{y) - figure 13- with a bias flux ^xstat also coupled to 
the ring. In previous papersiiSS this bias flux was used 
to control the behaviour of the SQUID ring alone or the 
ring interacting with one field mode. In the current pa- 
per it is used to control the interaction between two field 
modes via the non-linear properties of the SQUID ring. 
At first sight there might appear to be no a priori reason 
why, in this three mode system, it should prove easy to 
couple all the components together strongly. However, 
at least for the case of weak inductive coupling between 
the modes, we shall show that well characterized energy 
exchange can take place at (or close to) certain specific 
values of ^xstat- In these regions of bias flux multipho- 
ton absorption and emission processes occur. Thus, the 
energy required for an interaction to take place is approx- 
imately equal to the energy transfer in the absorption or 
emission of an integer number of photons with frequency 
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ojei in the first field mode to an integer number of pho- 
tons with frequency in the second. As for the two 
mode (ring+field mode) system, these are the exchange 
regions where the effective couphng becomes strong be- 
cause of the non-hnearity of the SQUID ring. As we shaU 
see, it is in these regions that many interesting quantum 
phenomena can be observed. 

To emphasize that the couphng across the extended 
three mode system is controlled by the bias flux, we cal- 
culate the average number of quanta in each mode and 
show that there is a large exchange of energy between 
the three modes at specific values of ^xstat, thus demon- 
strating that frequency conversion can take place in the 
system. In addition, we also calculate the second order 
correlation gl (again controlled by the bias flux) which 
quantifies the quantum statistics (bunching of quanta) 
for all three modest. We show that as the system evolves 
in time, strong entanglement occurs between the three 
modes. We quantify this by calculating various entropic 
quantities based on the von Neumann entropji2ii2Si2Si21. 
These are chosen for convenience and familiarity and be- 
cause they can be used to quantify the degree of entangle- 
ment between the subsystems (field modes and SQUID 
ring). Although these entropic quantities do have defi- 
ciencies as measures of entanglement^, there is no real 
consensus about which is the preferred measure within 
the quantum technology community. In the absence of 
any consensus, we opt for a familiar choice. 

II. THE TWO MODE HAMILTONIAN 

A. The SQUID ring in a classical field 

In our earlier work we treated the em field classi- 
calljiiiSS and the SQUID ring quantum mechanically, us- 
ing the well known Hamiltonian^ 

" ~ 2Cs 2As 

-fii/cos I 27r— I (1) 

Here, $s, the magnetic flux threading the ring, and Qs, 
the electric displacement flux between the electrodes of 
the weak link in the ring, are the conjugate variables^S*^ 
for the system (with [<I>s,Qs] = ih), fw/2 is the matrix 
element for Josephson pair tunnelling through the weak 
link and '^xem is the amplitude of the classical magnetic 
flux at the ring due to the em field mode. With ^xem set 
to zero we can solve the time independent Schrodinger 
equation to find the eigenvalues of the SQUID ring alone 
as a function of applied flux '^xstat/^o- As an exam- 
ple, we show in figure |21 the first three eigenenergies of 
the ring [£;„=o,i,2 {'^xstat/^o), where k = denotes the 
ground state, etc.] over the range < '^xstat/'^o < 1 
using parameters typical of a quantum regime SQUID 
ringii3£, i.e. = 1 x lO-^^F, A^ = 3 x 10~i°H 




External Flux $xski/*0 



FIG. 3: First three energy eigenvalues i?K=o,i,2 of a quantum 
mechanical SQUID ring as a function of bias flux '^xstat/^o 
over the range < ^xstat/^o < 1 for Cs = 1 x 1Q~^^F, 
As =3xlO-'^°H [hios = 0.043$o/As) and hu = 0.07<E>o/As = 
1.63hujs {u = 1.63lJs). 
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FIG. 4: First three time averaged Floquet energies as a 
function of ^xstat/^o for the SQUID ring of figure |3 where, 
again, TujJs = QM'i^l/Ks (ff ~ ITHz) and hv = 0.07$o/As 
[v = 1.63lJs) with a classical em field of frequencyuxem = 
O.Sojs and amplitude ^xem — 2 x 10~^$o applied. Here, the 
energy has been normalized to the ring oscillator energy hujs ■ 



(hence fiw, = 0.043$^/A, or |j = 9.8 x 10"Hz) and 

hv = 0.07$g/As = imnWs{v''=lMuJs). with <^xem 

turned on, we can use ^ to solve the corresponding 
TDSE. Again, by way of illustration, we show in fig- 
ure ^ the computed time averaged ring energy expec- 
tation values for the first three Floquet states (eigen- 
values of the evolution operator after one period of em 
field evolution) as a function of '^xstat/^o using the ring 
parameters of figure |31 Here, We has been set at 0.5ws 
with the associated ^xem = 2 x 10~^$o- As can be seen, 
energy exchange between these time averaged energies 
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occurs at specific values of the bias flux ^xstat and, as 
we have already pointed out, the number and position in 
^xstat of these exchange regions depends on uje, ^xem 
and the energy level structure of the ring. We have 
observed that these transition (exchange) points occur 
for values of bias flux such that (at least for small em 
field amplitudes ^xem.) NfnOf, « Ei{'^x) — Ej{^x) where 
N = 0,±1,±2, . . i We note that for the SQUID ring 
we can write down a renormalized oscillator frequency 
ris = LOs + 4:h'^n'^i^^Q^C~^ujJ^ which is related to the 
fact that there is an aja^ term in a Taylor expansion of 
the cosine (Josephson) term in the ring Hamiltonian(^. 



B. The SQUID ring in a non-classical field 

In the fully quantum description the Hamiltonian Ht 
for the SQUID ring-em oscillator mode system can be 
written as^*' 

Ht ~ H(. + Hs — Hes ■ (2) 

where and are, respectively, the Hamiltonian con- 
tributions for the field and the ring and H^s is the inter- 
action energy linking these together. 

Following equation the Hamiltonian for the 

SQUID ring alone is^ 

„ Ql , {^s - ^xstatf . 

^^=2C: + 2A^ ^'^^^^(.^-i^j- 

while the Hamiltonian for the em field (modelled as 
a parallel capacitance (Cg) - inductance (Ag) cavity 
mode equivalent circuit with infinite parallel resistance 

on resonance) takes the form Hf. = — h -^j^. Here, 
and Qe are, respectively, the cavity mode magnetic 
flux and charge operators for a field mode frequency 
We = 1/ \jC^Kf,. This cavity mode is coupled induc- 
tively to the SQUID ring with a coupling energy iJes = 
^{^s-^xstat)^e. where [i^^ is the em field-SQUID 
ring fiux linkage factor. 

By making a suitable transformation [using the unitary 
translation operator T = exp ^—i^xstatQslf^\ the Hamil- 
tonian ||3J) - now in calligraphic script - can be written 
more conveniently as 

(4) 

while the Hamiltonian for the em field mode remains un- 
affected. However, the interaction energy does transform 
to M'f^s — ^^s'^e- We denote the magnetic flux depen- 
dent eigenstates of by jcr)^. In our computations 
we then use a truncated energy eigenbasis both for the 
ring (|(T))and the em field mode (|n)). The basis states 
|(t) , |n), where cr = a, . . . , £7 and n = 0, . . . , A'^, are taken 
so that and A^ are, respectively, much greater than the 
average number of quanta in the ring and em field. 
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FIG. 5: Eigenenergies of the SQUID ring-em field mode sys- 
tem (normalized in units of hiDa) versus ^^atatj^o using the 
ring parameters of figure |21 with Cs = Ce = 1 x 10~^^F, 
clIs = 2lo^ (where ^ YTHz) and v = l.GSojj,. In (a) the 
field mode- SQUID ring flux linkage factor /lea = 0; in (b) 
/ies = 0.1. 



Using this truncated basis, we can then solve 



^ \in) = S„ |e„) (5) 

to obtain the eigenfunctions and eigenenergies of the two 
mode system Hamiltonian J^. The eigenenergies for the 
ring- field mode system are shown in figure for the C^, 
As and Kv values used for figure |21 with lOs — 2L0e, as in 
figure^ Here, in figure [SJa) the field mode-ring linkage 
factor /ies — 0.0 while in figure El^b) it is 0.1. From 
this spectral decomposition we can form the evolution 
operator via 



c/w = Ei^«)^^p(-^)(^"i (6) 

n ^ ^ 

The time averaged energy expectation values {{^i=s,e)) 
for the ring-*^ and the field (i.e. ,3^s and Mi) can then be 
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FIG. 6: Time averaged energy expectation values of Jifs 
(in units of hujs) as a function of ^xstat/^o for compari- 
son with figure 0] where the em-field mode is in a coherent 
state and the SQUID ring is in one of its first three energy 
eigenstates, i.e. the initial states are \a — i\/20)^ ® |o")5«, 
(cr = 0, 1, 2). Here, as for figure El = Ce = 1 x 10'^'^ F, 
LJe ~ 0.5tJs (where ^ ~ ITHz) , v = 1.63ljs but with fies = 
7.6 X 10~* for comparison with figure |1| 



calculated using the expression 



{{J^)) = lim - 



TY[p,{t)J^]dt 



(7) 



where pi, i — e, s are the reduced density operators for 
the em mode and the SQUID ring, respectively. In prac- 
tice, we have been able to ensure the convergence of 
the integral (O by integrating numerically from up to 
20,000/ws. 

Clearly, provided we select the correct initial state for 
the em field in the two mode system, there should exist 
a correspondence between the Floquet method used in 
section Hi Al [using |^] and the result of a fully quantum 
mechanical calculation. As an initial state the coherent 
state is an obvious choice since it is the closest quantum 
state to a monochromatic em field, as used in the Floquet 
approach above. With this choice we would expect a rea- 
sonable agreement between the fully quantum and quasi- 
classical computations. In fact the match between these 
two approaches can be very good. To compare with the 
quasi-classical result of section III Al we set the zero time 
(t = 0) product state as \a = i\/20)<f <8> \c^)s^, where \a)g 
is a coherent state of the field {ae\a)g = a\a)s )• Us- 
ing this coherent state we show in figure the calculated 
{{Jifs)) for an integration time t = 2 x IO'^/lJs with the 
energies normalized in units of fkUs- The computations 
have been made over the range < ^xstat /^o ^ ^ for the 
values of = 0, 1, 2, using the SQUID ring capacitance, 
inductance and $2;-dependent energy level structure of 
figure^. Here, as for figure 13 we have made uje = O.Sw^ 
while setting the flux linkage factor p^s = 0.00076. We 



chose this value of Hes so that the amplitude of oscil- 
lation of the coherent state in the em field coupled to 
the SQUID ring is equivalent to that used in the Flo- 
quet calculation of figure^ It is apparent that the time 
averaged energy expectation values, and their exchange 
regions, calculated in the quantum model as a function 
of ^xstat are very close to those found using the quasi- 
classical Floquet approach. To emphasise this, we show 
in figure 13 the two calculations superimposed. The very 
close fit between the two gives us confidence that the 
quantum model is physically valid with an accurate cor- 
respondence to the quasi-classical regime when a coher- 
ent state is chosen for the em field. 

As we have already noted for the SQUID ring in a 
classical field (section III A|l . we observed that transi- 
tion (or strong coupling) regions occur when NfiWe ~ 
Ei{^x) — Ej{^x)- From figure |31 we can see that this sit- 
uation exists when degeneracies in the spectrum of the 
system are lifted due to the coupling term between the 
Hamiltonians for the SQUID ring and em field. 



III. THE THREE MODE HAMILTONIAN 

Following on from Q , we now consider a SQUID ring 
with oscillator frequency Ws, threaded, as before, by a 
static bias flux ^xstat, but now coupled to two em field 
modes of frequency w ei and uje2 (see figure 0). With 
the usual flux ($i)-charge (Qi) commutation relation 
[$i,(5j] = «^<5y , the total Hamiltonian for this sys- 
tem can be written as 



(8) 



where the SQUID ring Hamiltonian is given by Q 
which has been transformed into the ^xstat basis. We 
choose to write our component Hamiltonians in ||SJ| in 

and cre- 



terms of the annihilation 



ation 



operators as 



tyc-e2 



hwe^ (flt^a 62 



where the position and momentum operators can be de- 
fined in terms of the magnetic flux and the charge op- 
erators via Xi — \J CiiOijh^i and "Pi = yJY/CihuJiQi for 
oscillator frequencies Ui = l/VCiAl, with the subscript 
i denoting 61,62 for the fields or s for the ring. Hence 
the Hamiltonians of the components of the system are 
identical (but extended to include an extra field mode) 
to those used in lllBI but written in terms of creation and 
annihilation operators. 
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FIG. 7: Comparison between time averaged energies of figures^ (red) and El (blue). 



The interaction energies J^e^s and J^se^ in ©, each 
of which represents the inductive coupHng between the 
SQUID ring and the oscillator modes ei and 62, respec- 
tively, are given by 




IV. TIME EVOLUTION OF THE THREE MODE 
SYSTEM 

The Hilbert space ^ for the SQUID ring-em field sys- 
tem is a tensor product of the Hilbert space .y for the 
ring and the Hilbert spaces S'l and S'2 for the fields, that 
is to say ^ — S\®S^(dS2. We denote in roman script the 
simple harmonic, em oscillator mode number eigenstates 
|n)i (a|^ai|n)i = n\n)i). In representing the SQUID ring 
we use greek script, for example |a)s, |/3)s, |7)s, • • ■ , to 
represent the eigenstates of the Hamiltonian in order 
to distinguish these from the number eigenstates \n)s of 
the ring {a\as\n)s = n\n)^. 

In dealing with the time evolution of the coupled three 
mode system we must first solve the eigenproblem [see 



©J . As in sectionnifor the two mode system, we use a 
truncated basis. This has the form 

{ |iV„,«,™) EE \n)g^ ® \k)^ ® \m)g^ I 
n = 0, ...,A^i, K = a, ...,r2, to = 0, ...,N2]. (9) 

where \n) and |to)^^ are the number states for the two 
field modes and are the energy eigenstates of the 

SQUID ring Hamiltonian. Here, iVi, 17 and N2 are taken 
to be much greater than the average number of quanta in 
each component of the system. With the eigenfunctions 
and eigenenergies of 10 , but using the three mode Hamil- 
tonian, the evolution operator can be calculated using the 
expression (|HJ. Then, assuming that the system at i = 
is described by the density matrix p(0), the density ma- 
trix p{t) = U{t)p{0)U''{t) at a later time t can be found, 
as can the reduced density matrices = Tr ^igj^a if) ^ 
Ps = Trsi(»S2 Pe2 = Tr^i®J^ (p) and pe,s = Tr^2 
Pse2 = Tr^i (p), Peie2 = Trj^(p). With these density 
matrices determined, we can then investigate a range of 
parameters which reveal much of the quantum behaviour 
of the three mode system, for example, the von Neumann 
entropy. 

In the following sections of the paper we present nu- 
merical results demonstrating various aspects of this be- 



haviour. Since the examples given are intended to be 
illustrative in nature, for simplicity we have made the 
capacitances for all three modes of the system the same, 
these being typical of quantum regime oscillators oper- 
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ating at a few K, i.e. Ce^ 
Amongst other things we wish to show that quantum fre- 
quency conversion can occur between the two oscillator 
field modes, via the SQUID ring. We have therefore made 
the two mode frequencies differ by a factor of two, i.e. 

= 2[Je2 while again, for simplicity, setting uj^.-^ = ujg. 
With all capacitances identical these frequencies corre- 
spond to — |Ae2 = Kg for a typical SQUID ring 
inductance A^ = 3 x 10^ ^"H (see section HB . Again, as 
in sectional we have put hv = O.OT^q/As {i^ = 1.63ujs) 
which is typical value of the pair tunnelling matrix el- 
ement for quantum regime SQUID rings operating at a 
few K. In addition, we have set the ring-field mode flux 
linkage factors at /Xe^s = 0.01 and ^se2 = 0.1 which ap- 
proximates to some reported experiments in the litera- 
ture involving two oscillator field modes coupled through 
a SQUID ring42-43. 



Strong Coupling of the SQUID ring to em field 
modes 



In sectionniwe demonstrated that strong coupling be- 
tween the SQUID ring and the em field occurs when 
degeneracies in the spectrum of the Hamiltonian are 
lifted due to the coupling between the components of 
the system. Numerically, this equates to the condition 
Nhuje ~ Ei{^x) — Ej{^j;), where N is integer and the 
Ei{^x) and Ej{^x) are the and j*'' eigenvalues of the 
system Hamiltonian Jift- These regions of strong coupling 
(the exchange regions), which develop at specific values 
of the bias flux {^xstat) applied to the SQUID ring, are 
dependent on the ring eigenenergy structure and the fre- 
quency of the em field. Similarly, for the three mode 
system strong energy exchange will occur between the 
two field modes and the SQUID ring when the coupling 
terms lift degeneracies in the spectrum of the Hamilto- 
nian. This will occur when 



^,($,) -£,($,) (10) 



where, now, the Eij{^x) are the i'^'^/j*-"- eigenvalues of 
J^y and -/Vei,e2 are integers (positive or negative). In 
this case N^.-^ photons with frequency ujei in the first 
field mode are used to excite the SQUID from the j*^- 
state to the i'^'-state with the emission of N^^ photons 
of frequency uJe2 into the second field mode. Taking the 
ring-field parameters set above, we have calculated the 
energy eigenvalues for the three mode system. These are 
shown in figure|Hl It can be seen that these eigenenergies 
possess a very rich structure which leads directly to the 
results presented in this paper. 

Given a choice of the initial state for the fields in the 
two oscillator modes, we can calculate the time aver- 
aged energy expectation values (normalized in units of 
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FIG. 8: Eigenenergies (in units of hus) versus ^xstat/^o 
of the three mode (em mode-ring-em mode) system with 
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tors jj,ei3 and fj,s 
while in (b) fj,eit 



C. 



= 1 X 10" 
ITHz), v — l.GSujs and flux linkage fac- 
. In (a) /ieis and /^se2 are set equal to zero 
= 0.01 and /isej = 0.1. 



hws = h/ VCsAs) of J^ei , and Jifg using equation Q 
and integrating numerically from up to 2 x lO^/wg. As 
an example, we show in figure |5fa) these time averaged 
energies plotted as a function of ^xstat/^o assuming that 



ait — Q the system is in the state |1) 



®|0), 



I.e. 



with one photon in the first em mode and none in the sec- 
ond. The peaks are a manifestation of the strong coupling 
between the various oscillators at and around specific val- 
ues of ^xstat I 'I'o where the energy transfer between the 
various components of the three mode system occurs. 
Thus, starting in the initial state |1)^^ ® W) ^ ® lO)^^, 
it can be seen that in the exchange regions, on average, 
energy is being transferred from the first mode to the 
SQUID ring and to the second mode. This will become 
more transparent when we compute the time evolution 
of the expectation values of the number operators for 
the components of the system (below). In figure Efb) we 
show one of the three mode exchange regions (around 
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FIG. 9; (a) Time averaged energy expectation values (in units 
of huos) over the range < ^xstat/^o < 1 for the three mode 
system of figure |S| and (b) as for (a), but expanded around 
^xstat/^o = 0.426. The initial state for the system is |1)^ ® 
® \Q) g and the circuit parameters are as in figure |H| 



^xstat/'^o = 0.426) of figure a) on a much expanded 
scale so that the details can be seen more clearly. 



B. Quantum statistics of the SQUID ring-field 

system 

An important aspect of non-classical electromagnetic 
fields is the quantum statistics of photons (bunching of 
photons) described by the second order correlations^SiS 



9i 



(A^f)-(A^,>, 



(iVf) =Tr[p, (ala,) 



Ml 



ei, 62, s 



(11) 



The value of g^^^ = 1 corresponds to Poissonian statis- 
tics. Values of g^^^ greater than 1 indicate photon bunch- 
ing (i.e. where the photons arrive in groups) while values 
of 5^^'' smaller than one indicate antibunching (i.e. the 
regular arrival of photons). The latter regime is charac- 
teristic of non-classical electromagnetic fields since it can 



be shown that in classical optics (7*^^^ > 1. Only in quan- 
tum mechanical systems can g'^'^'^ < iMiM. This is well 
known in quantum optics but is, perhaps, less familiar in 
condensed matter physics. 

In this work we show that the statistics of photons 
threading the SQUID ring affects the statistics of elec- 
tron pair condensate tunnelling through the Josephson 
junction in the ring. This is quantified by the second 
order correlations although higher order correlations can 
also be calculated22i. These correlations fully describe the 
quantum statistics and quantum noise of the photons in 
the two field modes and the superconducting condensate. 



C. Quantum Entanglement in the Three Mode 
System 

The creation of entangled states of multi-particle sys- 
tems is a key feature of all quantum technologies. In their 
pursuit the generation of entanglements in real physi- 
cal systems is clearly of very considerable interest. In 
this regard, it appears that the non-linear properties of 
the SQUID ring can be used very efhciently to entan- 
gle circuit subsystems (here, field oscillator modes) that 
are coupled to it. As we shall also show, the ring non- 
linearity can also be used with facility to generate energy 
conversion between the two oscillator field modes. Again, 
taking as our example the three mode system of figure |21 
we shall demonstrate that as this system evolves in time 
its three components become, to a greater or lesser ex- 
tent, entangled. The degree of this entanglement can 
be quantified by using entropic quantities. The entan- 
glement for a two mode (ring-oscillator) system can be 
quantified bji^lS*^ 

Tab = S (pa) + S (pb) - S (p) (12) 

where S {p) is the von Neumann entropy given by 

5(p) = -Tr[pln(p)] (13) 

with Pa — TtbP and ps — Tr^p. This entanglement en- 
tropy is positive or zero (subadditivity property of the en- 
tropy) . Examples of the calculation of this entanglement 
for the two mode system can be found in our previous 
workSS. 

An analogous quantity can be used to characterize 
the entanglement of a three component system^SiS^i^^iS. 
Thus, for the field mode-SQUID ring-field mode system, 
this takes the form 

I = S (p.^, ) + S (ps) + S (p ) - S (p) (14) 

which can be written as 

I = leis + Ise2 + I ieis;se2) (15) 

where /g^s and /se^ are the entanglement entropies be- 
tween (01 (g) and y £"2, as defined above l(T^ . and 

/(eis; 862) = S (pets) + S {pse2 ) - S [ps) - S (p) (16) 



9 



describes a deeper entanglement between S'l (g) y and 
(g) (o2. Understanding of this deeper entanglement, 
that exists in three component systems, is intimately 
connected to the strong subadditivity property of the en- 
tropy. This can be used to demonstrate that the quantity 
/(eis; 562) is positive or zero. We note that to prove this 
presented a very difficult problem in the theory of en- 
tropy (it was a conjecture for many years until a proof 
was providedSi^i^iiS). 

In this paper we are not going to proceed further into 
this deep problem of entanglement in three component 
systems as it has been discussed in detail elsewhere^. 
However, we note that the entanglement / can also be 
expressed as: 



/ =5(peJ+5(p.)+5(peJ-5(p) 
= leis + Ieie2 + I (giS; 6162) 
= Ie2S + Iete2 + I (e2S; 6162) 



(17) 



where /(eis;eie2) and /(e2s;eie2) are non-negative 
numbers. On physical grounds, i.e. because we are us- 
ing the SQUID ring as the intermediary between the two 
field modes, we choose to show numerical results for the 
entanglement between the SQUID and the first mode 
{leis), the SQUID and the second mode {le^s) and also 
the /(eis; 562). 
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V. NUMERICAL CALCULATIONS 

As we have shown in figures El and ^ strong coupling 
between the various of components of ring-field mode 
systems only occurs over small regions in ^xstat - the 
exchange regions. We now see how the variation in cou- 
pling across an exchange region affects the number oper- 
ator expectation values, quantum statistics and entangle- 
ments - all important quantities reflecting on the quan- 
tum behaviour of these systems. Continuing from fig- 
ure El we calculate these quantities at each of the three 
flux bias points A, B and C (at ^xstat/^o = 0.0,0.4246 
and 0.4263, respectively). In each of the following com- 
puted examples we assume that at t = the first field 
mode (ai contains one or more photons while the second 
contains none. In our first set of examples we choose the 
t = state in the first mode to be a number state; in the 
second set we make this a coherent state where 
ctei\^)si = For the case of the number state 

we assume that at t = the three mode system is in 
the state ^ ® |a) ,y ® |0)^2 ■ For the example where we 
adopt a coherent state for the first mode we choose for il- 
lustrative purposes (and computational ease) the system 
state \A = iy/3)si ® 1").^ ® lO)^^- 

A. Number state computations 

As is evident from figures OJa) and Elb), the flux 
bias points have been selected either to be well away 



FIG. 10: Starting in a pure state |l)<fj^ (g)|a)^®|0)^2 5 plots of 
(a) the number expectation values (riei), (ws), (^62) and (b) 
the entanglement entropies /eis, /se2) /(eis;se2), / versus 
dimensionless time uust for the three component system of 
figure |5| with the ring flux biased at A in flgure|5Ia). Here, 
the system parameter values are as for figure |HJb). 



from, or within, an exchange region, i.e. point A and 
points B and C, respectively. For a complete, quantita- 
tive view of the system we should compute, in sequence, 
the number expectation values (nei, 5,62)7 the entropies 

{leis 1 Ise2T I i^i^'j ^^2) and /)and the (7ei's,e2 correlations 
for the first held mode {S'l), the S QUID ' ring (^) and 
the second field mode {(02) as a function of normalized 
time LUgt. However, it is apparent in figure [TUI (bias point 
A) that, starting in a pure state (S> \a}y d) \0)s2 

for S'l, the number expectation values (a) and entropies 
(b) remain constant as a function of time. We note that 
(us) is not zero because the ground state \a)y, of the 
SQUID ring is not the same as the ground state of a sim- 
ple harmonic oscillator |0). From the definition given in 
section HV Bl the fact that (ngj) — makes the calcula- 
tion of g*^^^ for the second mode at bias point A physically 
unmcaningful since this is division by zero. It is there- 
fore very sensitive to numerical error. However, such is 
not the case when the bias point point is shifted into an 
exchange region. Starting again with the system state 
|l)(?i (8) |a)^ ® |0)^27 we show in figures [TTlfbias point B) 
and 1 1 21 f bias point C) the average number of quanta - the 
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FIG. 11: With the initial pure state |l),fi (8) |a) <8) |0)^2 , and 
with the parameter values of figure |HJb), a plot of the num- 
ber expectation values {riei), {»^s), (fie2) versus dimensionless 
time ujst for the static magnetic flux on the SQUID ring set 
at point B in figure l^^b). 



FIG. 13: Entanglement entropies leis, Ise^, I (eis; 562) and I 
versus dimensionless time uigt for the three component system 
of figure]^ with the parameter values of figure |HIb), starting 
in state \l)ei \a) <S) \0) ^2 with the SQUID ring flux set at 
point B [figure Elb)]. 




1000 2000 3000 4000 

Dimensionless Time (w^t) 



5000 



FIG. 12: Starting in the pure state ® |a) ^ (8) |0)^2i 
with the parameter values of figure |SJb), plots of the num- 
ber expectation values (riei), ('^s), (j^ea) versus dimensionless 
time uist for the three component system of figure |n| with the 
SQUID ring flux biased at point C in flgure|UJb). 



(rii) - in each of the three modes as a function of time. 
With this choice of starting state these results demon- 
strate clearly the quasi-periodic exchange of energy be- 
tween the various components of the system. Since the 
exchange coupling is strongest at C, this is where we 
would expect to find the maximum energy transfer be- 
tween the first and second field modes, as is the case 
(figure El- We note that in the computed results of fig- 
ure 1121 the second field mode number expectation value 
(and that of the SQUID ring to a much smaller extent) 
is a maximum when that for the first field mode is a 
minimum. This is the signature for frequency (down) 
conversion, in this example by a factor of two. The pro- 
cess could, of course, be run backwards to generate fre- 



quency up conversion from the second to the first field 
mode via the quantum non-linearity of the SQUID ring. 
Given that this non-linearity can be to all orders, we see 
no obvious reason why much higher ratio frequency con- 
versions should not prove practicable. 

From a theoretical viewpoint the problem with demon- 
strating high ratio frequency conversions is the rapid rise 
in the number of basis states required as the down (up) 
conversion frequency ratio increases. The computational 
difficulties increase accordingly. Nevertheless, even given 
the limitations on the computational power we have avail- 
able (Compaq XPIOOO alphaserver with 2GB RAM), we 
have been able to demonstrate quantum down conver- 
sion by a factor of ten in frequency. We intend to deal 
with this in a future publication. There is some indi- 
cation that these down conversion processes occur—. In 
figurelT^the input state is the number state \ but, as 
we shall show, down conversion can occur for a coherent 
input state. It may well be that this ability to generate 
photon down/up conversion could have practical applica- 
tion for pure state sources in quantum information pro- 
cessing and quantum computing. For example, it may 
prove desirable to take single photon terahertz sources, 
as are now being developed'**'*^, and use these to provide 
the input state to a SQUID ring to generate photons at 
much lower frequencies suitable for solid state quantum 
circuit technologies. It is also clear that if very large 
down/up frequency conversion ratios can be achieved ex- 
perimentally, there could well be interesting metrological 
applications, for example, in frequency standards. 

This role in linking the two field modes together in a 
strongly non-linear, quantum mechanical manner is em- 
phasized in figures ^1 and ^1 . Here, the time vary- 
ing entanglement entropies are computed for bias points 
B (figure [T5|) and C (figure following the definitions 
given in section llV CI Here, again, we have started the 
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FIG. 14: Entanglement entropies leis, Ise2, I {^is; se2) and / 
versus dimensionless time ujst for the three component system 
of figure |5| with the parameter values of figure |HIb), starting 
in state \ (S)\a) ^ (gi\0) with the SQUID ring flux biased 
at point C in figure I^Jb). 



system in state (^\a),y |0)^2- I* can be seen 

that the entanglement between the various components 
of the system (leis, Ise2 and /(eis;se2)), and the total 
entanglement entropy for the system (/), are stronger at 
C than at B which, from figures l^b). 1111 and 1121 is to 
be expected. In our opinion it is this ability to control 
the degree of entanglement between the components of 
(for example) this three mode system simply by changing 
^xstat which marks out the SQUID ring as a potentially 
very useful device in future quantum circuit technologies. 
This is emphasized by the contrast between figures lTUT b') 
(for bias point A) and where the system and sub- 
system entanglement go from zero to a maximum for an 
adjustment in ^xstat around 0.04<I>o. 

Although, in the above we have considered in some 
detail entanglement between two field modes (input and 
output) interacting via a SQUID ring, there are many 
other coupled systems of field modes and SQUID rings 
which could be studied. One which may be of impor- 
tance, both scientifically and technologically, is an input 
field mode linked through a SQUID ring to two sepa- 
rate output modes at half the input frequency. From 
the results obtained in this paper we would expect the 
two (down converted) photons to be strongly entangled 
with the degree of entanglement controlled, again, by the 
bias flux ^xstat applied to the ring. Furthermore, given 
the non-linearity of the SQUID ring, we would also ex- 
pect it to be possible to entangle a large number of out- 
put photons starting at an initial input frequency and 
down converting to a whole set of lower frequency output 
modes. As a technique, the use of a SQUID ring to gener- 
ate entanglements between several systems could well be 
applied to great advantage in fundamental experimen- 
tal studies of quantum mechanica^SiiSiS. It could also 
have implications in quantum computing, for example. 



in creating an entangled input register for a quantum 
computer. It has also been suggested that the creation 
of (large number) multi-particle entangled systems could 
lead to new sensors and instrumentation of unparalleled 
sensitivity'*'' and it may be that SQUID rings are very 
well suited to creating these entanglements, at least for 
photons. 

There are other possible ways that the input register 
of a quantum computer could be based on the non-linear 
properties of SQUID rings described in this paper. For 
example, we could set the three modes of the coupled 
system in figure [21 all to have the same oscillator fre- 
quency Us- Then, with ^xstat biased within an exchange 
region, we could arrange to create a qubit superposition 
state of |0) and |1) in the output mode starting from 
the number state |1) of the input mode. As our results 
have demonstrated, this could be done in such a way 
as to ensure that the input and output oscillator modes 
are entangled. Once the desired qubit state of the out- 
put mode had been realised, in principle the bias flux 
could then be switched away rapidly from the exchange 
region (or switched off) thus leaving the input and out- 
put modes entangled but uncoupled. An array of these 
circuits could then be used as an qubit register for a quan- 
tum computer, where the qubits would be entangled but 
not coupled to the input modes. Conceivably, this ar- 
rangement could facilitate quantum error correction for 
quantum computation*®. Schemes of this kind may well 
find application in quantum encryption and transmission 
of information at a more complex level than is usually 
considcrcd^i^Si^. 

Since the underlying purpose of this work is to demon- 
strate the influence of the SQUID ring non-linearity on a 
coupled quantum system, it is important to show quan- 
titatively the way in which the quantum statistics of 
the photons affects the quantum statistics of the elec- 
tron pairs (i.e. the superconducting condensate flowing 
through the weak link in the ring). As we explained above 
f section llVB|) . this is quantified with the second order 

correlations g^^'^K In figures [T^ and ITCl we plot the second 

order correlations gi'f^ , gi^^ and gi'^'^ for bias points B and 
C as functions of time with, again, a starting state for the 
system of |l)^j (^\a),y \^)s2- "with the number ex- 
pectation values and the entanglement entropies, we see 
a strong oscillatory behaviour, particularly in figure 1161 
(biased at point C). In order to interpret the results we 
first note that for the state \Q)i {i ^ S'l, (02, S^), or other 
states close to this state, the average number of photons 
is near zero. From equation Hll|) this means that the sec- 
ond order correlation becomes very large. We also note 
for the number state |1), {N'^) — {N) = and the corre- 
sponding g^^'> — 0. With this in mind, we see in figure [T31 
(bias point B) that the first field mode, in number state 
|l)^j at i = 0, starts with a g^^^ at zero. It remains 
extremely close to this value over the time of the compu- 
tation, i.e. this field mode stays reasonably close to the 
number state jl),?!. By contrast, the second field mode. 
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FIG. 15: Second order correlations g^s^ ,g^^^ and g^^^ versus 
dimensionless ujst for the three component system of figure|^ 
starting in a pure state \l)gi ® \ct)^ <8> |0)^2 with the static 
magnetic flux on the SQUID ring set at bias point B in flg- 
urel^Jb). Here, the system parameter values are as for flg- 
urelHIb). 
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FIG. 17: Starting the three component system with the first 
field mode in a coherent state, i.e. \A = ix/S),?! ® \q) ^ ® 
1 0)^2 1 ^iid using the parameter values of figure IHTb). plots 
of the number expectation values (riei), (ws), {ne2) against 
dimensionless time uJat for the static flux on the SQUID ring 
set at point C in figure I^Jb). 
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FIG. 16: Second order correlations (?ei , 5s and ge2 versus 
dimensionless uigt for the three component system of figure |S| 
[parameter values as for figure |HJb)], starting in a pure state 
1 1)^1 ® ® |0),g'2 with the SQUID ring bias flux set at bias 
point C in figure |UJb). 



assumed to be in the state |0)^2 at i = 0, has on average 
a large g^"^^ value, although this dips well below unity 
almost periodically with time. This demonstrates that 
even at bias point B, on the edge of the exchange region, 

(2) 

the value of regularly falls below one. It also shows 
that for this initial condition the quantum statistics of 
the second field mode cannot be described by classical 
means. In figure ITCl (bias point C), at i = 0, we again 
assume that the first field mode is in state with 
the second in state \0)£'2- As before, the first field mode 

(2) 

starts at ge^ = but as the wavefunction for the system 

(2) 

evolves with time we see that gei regularly shifts away 
from zero. Correspondingly, the first field mode is no 
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FIG. 18: With initial state \A = iV3)s-, \a)s^ ® \0)g2 for 
the three component system, the entanglement entropies Zeis, 
Is 32, I (cis; 862) and I versus dimensionless time uist, where 
the system parameters are as for figure |HIb) and the SQUID 
ring is fiux biased at point C in figure |5Jb). 



longer in the pure number state due to its interac- 
tion with the rest of the system. The first field mode is 
described by a reduced density operator which, at these 
points, represents statistical mixture of states with a low 
photon number expectation value (as can be seen from 

figure [T^. As a consequence, f/e?^ increases since the de- 
nominator in equation Ijlll) becomes very small around 
these points. 



B. Coherent state computations 

In figures [T7I ll8land[T^we show the number expecta- 
tion values, the entanglement entropies and the g*-^-* cor- 
relations for the bias point C in figure |3b), taking the 
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FIG. 19: Second order correlations g^-^ ,gs and g^^ ver- 
sus dimensionless time uist for the three component sys- 
tem [parameter values as for figure |SIb)] starting in state 
\A = iV3) ,?i (g) |a)j^ (g) |0)^2 with the flux on the SQUID 
ring set at bias point C in figure Elb). 



initial state of the system as \ A — i\/3)s>^ (g) \a) ^<E'|0)^2, 
i.e. where the first field mode is in the coherent state 
\'iV3)si at t = . It is apparent (figure [T7|l that there 
is energy transfer, via the SQUID ring, between the first 
and second field modes of the system, just as in figure [T^ 
for the pure state d) \ct)y \0)s2 ■ Thus, when 

(riei) decreases (?T.e2) increases, and vice versa. How- 
ever, it is evident that the regular oscillatory behaviour 
seen in figures ITTI and IT^ has been lost. We note that, 
as in figures 1111 and 1121 the number expectation value 
of the SQUID ring remains at a roughly constant value, 
highlighting the view that the SQUID ring is acting as a 
non-linear control medium linking the two quantum field 
modes together. In figure^lwe see that the components 
of the system again entangled very strongly but, unlike 
the previous computations of figures ^1 and ^1 there is 
no longer any quasi-periodic disentanglement to be seen. 
The system remains entangled at all points in time, i.e. 
the total entanglement entropy is always high. In fig- 
ure 1191 there is clearly a significant deviation from the 
behaviour of the g^^^ coefficients displayed in figures 
and 1161 Thus, in figure [T^ is close to one for all of 
the time evolution (and for most of the time just greater 

than one) whereas gs spends most of its time just less 

(2) 

than one and ge2 is almost always greater than one. 

In all the examples given above we observe that, as 
expected, the system displays strong entanglement when 
the first field mode has a low number expectation value, 
i.e. the components of the system have evolved away from 
their initial pure states. We also note that when we start 
the first field mode in the pure state \l)g^(^\a),y(^\0)s2, 
at bias points B and C, we find strong entanglement 
between the various components of the coupled system 
which, at certain times (semi-periodically), disentangle. 
The entanglement entropies are, of course, theoretical 



quantities that demonstrate the development of quantum 
correlations between the three modes. In principle, ex- 
perimental observation of these quantum correlations can 
be achieved by determining Bell type inequalities^SiiSi^. 
In the context of the work presented here this will require 
further theoretical investigation. In figures 1 1 til and 1191 
it is evident that the g*^^' correlation coefficient for at 
least one of the components of the system becomes less 
than one at some point in the evolution of the system. 
Hence, for the initial conditions used in this work, we 
conclude that the photon statistics of this system cannot 
be described by classical optics. 



VI. CONCLUSIONS 

In this paper we have studied the coupling of a SQUID 
ring to two em field modes. For this we have made the 
assumption that the ambient temperature of the system 
is low enough to be able to treat each of the three com- 
ponents (ring -|- two field modes) quantum mechanically. 
Our purpose has been to demonstrate that the SQUID 
ring, as a non-linear quantum object, can be used to cou- 
ple a number of quantum oscillators together to generate 
physical phenomena of great interest. As we have em- 
phasized, there exist obvious parallels with the field of 
quantum optics. However, in quantum optical systems 
the coupling media involved generally display only weak 
polynomial non-linearity. In contrast the SQUID ring, 
with the cosine term in its Hamiltonian description, can 
generate non-linear interactions to all orders. It is this, 
plus the <I>o-periodic nature of its behaviour as a function 
of external flux which makes the SQUID ring of such 
interest in the burgeoning field of quantum circuit tech- 
nology. Viewed from the perspective of quantum optics, 
the SQUID ring (or a set of coupled SQUID rings) can be 
thought of as a non-linear medium par excellence which 
can easily create very strongly coupled regimes (albeit at 
lower frequencies - for example, at THz frequencies and 
below) which are inaccessible using conventional optical 
materials. 

In our theoretical investigations of the two em field 
modes coupled through a quantum SQUID ring we have 
also applied a static external magnetic flux {^xstat) to 
the ring. In this paper we started with the simpler exam- 
ple of a SQUID ring interacting with a single em mode. 
We showed that the coupling between the components of 
the system can be strong. This strong coupling only oc- 
curs over small ranges in ^xstat-, centred around speciflc 
values of this bias flux, i.e. in what we term exchange 
regions which are govern by the energy eigenstructure of 
the system. It is in and around the exchange regions 
that the quantum non-linear nature of the SQUID ring 
is made manifest through the coupling of the field modes 
via the ring. For example, in these regions energy can 
be exchanged between the field modes and the SQUID 
ring and, through the intermediary of the ring, between 
the field modes themselves. From this we suggest that. 
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suitably driven, this system may act as a frequency con- 
verter suitable for operation up to the THz range. From 
our viewpoint this illustrates the utility of these exchange 
regions since it is here that the strong quantum couplings 
develop between the components of the system. We have 
added to this perspective by calculating other physical 
phenomena associated with such a coupled quantum sys- 
tem. To illustrate this we have computed the statistics 
of the various quanta (quantified with the second order 
correlations) and the degree of entanglement (quantified 
with various entropies) between the components of the 
system, both outside and within the exchange regions. 
Our results demonstrate quite clearly that the mesoscopic 
SQUID ring can be used as a flux tunable element to ma- 
nipulate these (and presumably other) quantum proper- 
ties of these coupled circuit systems. As such, the work 
presented here may be of considerable relevance to cur- 
rent experiments on quantum superposition of states in 
SQUID rings and on probing crossing/ anticrossing re- 
gions of their energy level structureiSiiSiS. 

We note that we have neglected dissipation in the re- 
sults presented in this paper. Thus, we have computed 
the time evolution of the three mode system using the 
equation dtp = —i[H,p\. A more realistic calculation, 
with dissipation due to the environment taken into ac- 
count (for example, in the references 'SS^'Sf^/BCf) is now 
being developed. The equation for the time evolution 
then takes the form dtp = —i[H, p] + f{p), where a dis- 



sipative term f(p) has been introduced to represent the 
environment. We intend to extend our work to investi- 
gate in detail the effects of this environmental dissipation 
on the behaviour of SQUID ring-field mode systems. 

We suggest that the three component system (SQUID 
ring + two field modes), and its extensions, is rich in pos- 
sibilities for device applications (e.g. in quantum gates, 
quantum encryption and frequency conversion) ; it is also 
a pointer to more sophisticated quantum technologies in 
the futuraSi. Given that the technical problems associ- 
ated with such technologies can be overcome, it seems 
likely that the SQUID ring (and related weak link cir- 
cuits) will, in the future, be able to operate at THz fre- 
quencies. This could complement the current drive to 
develop THz applications in (classical) communications 
and imagin g^*^!^^ . It would also allow for quantum circuit 
technologies to be utilized at quite accessible tempera- 
tures. 
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